Radiative corrections to the three-body region of the Dalitz plot of baryon 
semileptonic decays with angular correlation between polarized emitted 
baryons and charged leptons: The initial-baryon rest frame case 

C. Juarez-Leon/ A. Martinez,' M. Neri/ J. J. Torres,^ Ruben Flores-Mendieta,^ and A. Garcia"^ 

^ Escuela Superior de Fisica y Matemdticas del IPN, 
Apartado Postal 75-702, Mexico, D.F. 07738, Mexico 

^ Escuela Superior de Computo del IPN, 
Apartado Postal 75-702, Mexico, D.F. 07738, Mexico 
^ Instituto de Fisica, Universidad Autonoma de San Luis Potosi, 
Alvaro Obregon 64, Zona Centro, San Luis Potosi, S.L.P. 78000, Mexico 
Departamento de Fisica, Centro de Investigacidn y de Estudios Avanzados del IPN, 
Apartado Postal 14-740, Mexico, D.F. 07000, Mexico 
(Dated: November 27, 2008) 

Abstract 

We complement the results for the radiative corrections to the S2 • 1 angular correlation of baryon 
semileptonic decays of Ref. Ij^ with the final results in the rest frame of the decaying baryon. 
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In a recent paper [1] we obtained the radiative corrections (RC) to the Dalitz plot of the 
semileptonic decay of a spin- 1/2 baryon A, 



Ao{p,) ^ B+ip2) + £-(1) +V,ip,), 



(1) 



when the angular correlation §2 • 1 between the spin §2 of the emitted baryon B and the direction 1 
of the emitted charged lepton i is observed. As is customary, the results were presented in the rest 
frame of B where p2 = (M2, 0, 0, 0). However, due to experimental conditions, it may be more 
convenient to produce such RC in the rest frame of A where pi = (Mi, 0, 0, 0). It is not possible 
to translate directly the final result of Ref. IIjOQ into the final result of the latter RC. The calculation 
must be retaken starting at earlier stages. In this paper we shall complement the analysis of Ref. ill] 
and present the final result for the RC to the §2 • 1 angular correlation in the rest frame of A. 

We shall follow the same procedure and use the same conventions and notation of Ref. 111]. We 
shall omit a detailed discussion, which can be found in this reference. Let us just recall that our 
results for the virtual part will be model-independent, gauge invariant, finite in the ultraviolet, and 
will contain the infrared divergence. In order to avoid repetition of long expressions, it will be 
convenient to trace a close parallelism with the analysis of the RC for the angular correlation Si ■ 1 
between the spin Si of A and the direction 1 of £ which were obtained in Ref. 121]. 

Without further ado, the differential decay rate with virtual RC in the rest frame of A including 
the §2 • I correlation and covering the two charge assignments of is 



dT 



V 



dT'y - dVP. 



The unpolarized part was already calculated in Ref. 121]. It has the form 



dT' = dVL 



a 



(2) 



(3) 



The full expressions for Aq, A'^, A'(, (p, and (p' can be found in Eqs. (Bl), (B2), (B3), (11), and 
(12) of this reference. Here the phase space factor dn = {l/2){Gl/2)dE2dEdned(j)22Mi/{27if 
is one-half the phase space factor of Eq. (25) of this same reference. The polarization appears in 



dT 



dQ 



a 



+ - {B'^ 

71 



'B'- 



S2-1. 



(4) 



§2 is introduced with the spin projector S(s2) = (1 — 75 ^2)/2 applied to the ub spinor. The four 
vector S2 obeys S2 ■ S2 = —1 and S2-p2 = 0. The trace calculation will lead to products S2 ■ a with 
a = l,pu,pi. These products are specialized to the rest frame of A using the relation JSI] 



S2 ■ a 



'P2 




'P2 / 


— ao - a 

.e/2 


= sf . 


M2 I 



ao - 



P2 • a 



— a 



(5) 
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which corresponds to the Lorentz transformation from the rest frame of B to the rest frame of A. In 
the first equality of ([5]) it is understood that p2, s^, and the components of the 4- vector a = (ao, a) 
are specialized in the rest frame of A. In the second equality of ^ only sf still remains in the 
rest frame of B. The upper indices A and B on the spin of the emitted baryon emphasize this fact. 
In addition, the contribution of the correlations sf • f>2 and sf • 6,^ to the sf • 1 correlation over 
the Dalitz plot are taken into account with the substitution rule 121] • p ^ (sf • i)(p • 1), with 
p = P2, Py, which is valid under integration of the variables, other than E and E2, contained in 
dQ. These steps will be extended to the bremsstrahlung part, where p = k will also appear. 
The results in Eq. Q are new. Their explicit expressions are 



A 



(^) 



Q'l 



Ml 
Mo 



+ Q2 



P22/0 



E 



B' 



M2 V E2 + M2 

-D,E'j + DsEip2yo + l), 
B" = -D3EI ■ p., 



(6) 

(7) 
(8) 



where 



Q'l = (F^ - Gl 



EE2 - P2 • 1 



Mo^ 



+ {F1G2 + F2G1 
— {F1G2 — F2G1 



Ml Ml 
EM^ + {E-El){EE2-V2-\) 



Ml 

M2{EE2-p2-i + m^] 



Mi 



m 



+ (FiGa + FsGi)j^iEE2 - P2 • 1 + M^) 



— {F1G3 — F3G1 
-F2G 



+ (FiF2 + GiG2)- ' ^ 



Ml ' ^ ' ' ' ' " Ml 
2^(-/^ - P2 • 1 + EEl) - (F2G3 + F3G2) 



— F'^G; 



'Ml 
M2m2 



Mf 



(EE° + /' + P2-l), 



(9) 



Q2 



iFl^Gl)^^iMi-E2)-iFl-Gl)^^'''- 
E~ 



'- + 2FiGi^^{El-E) 



+ (F1G2 + F2Gi)-^{El -E)- (F1G2 - F2Gi)-^(E;; - E) 



{F1G3 + F^Gi 



Ml 
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(F1G3 - F^Gi)^^ - (F1F2 + GiG2)^, (10) 
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D-^ = 2{f[g[ -g'l^), and D4 = 2{f[g[ + g['^). In Eqs. and hereafter and / will denote the 

magnitudes of the corresponding 3-momenta. To avoid making the notation more cumbersome, 
we did not put primes on the form factors on the r.h.s. of Eqs. ^ and (fTOl) . However, it must be 
kept in mind that it is the primed form factors f[ and g[, where all the model-dependence has been 
absorbed, that appear in these equations. We have limited ourselves to put primes on Q[ and Q'2, 
as a reminder of this fact. 

To the virtual RC of Eq. ([3]) one must add the bremsstrahlung RC. It arises from the radiative 
decay 

A-^ B + £ + iye + -f, (11) 

where 7 is a real photon with 4-momentum k = (u;, k) and which, in order to regulate the infrared 
divergence, is emitted with mass A and with an additional longitudinal degree of freedom. The 
summation over its polarization is performed according to Ref. [4] and its model-independent 
contribution is controlled with the Low theorem [|50. The integrations over k are performed 
covariantly following Ref. [4]. 

Introducing the projector S(s2) in the radiative decay transition amplitude, following the usual 
steps of squaring it and summing over all polarizations including the summation over the 7 
polarization j^, and extracting the ■ 1 correlation as explained in the virtual part, one obtains 
for the differential decay rate the result 



dr 



B 



dr'n - dT 



B ' 



(12) 



where dT'^ is independent of sf and can be identified with one-half the unpolarized decay rate of 
Eq. (27) of Ref. [|20. The B spin-dependent part is given by 



rfrg^ = -dn {loB' + Dsps + /^4P4) • 1. 



a 



TT 



(13) 



Io{E, E2) is the infrared-divergent integral given in Eq. (26) of Ref. [|20, B' is identified with 
Eq. dV]) of this work. The contributions which are different with respect to the corresponding ones 
of Ref. m are 



P3 



— [ dx 

^ 1-I3x 
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(14) 
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and 



P4 



P2[ 



dx 



dy 



2lT 



D{l-I3x) 



l-(3x 



l-Px 



- 1 



UJ 

E 



X 



where 



and 



D = E 



UJ 



^ + (p2 + l)-k, 

P2l{yo - y) 

D 



(15) 



(16) 



(17) 



The neutrino energy is E^ = E^-uj, where E^ = M1-E2- E, and yo = {E^ -pl-f) /2p2l. 
In Eqs. (fT4l) and (fT5l) k is refered to a coordinate axis system where 1 points along the 2;-direction 
and p2 lies on the (x, z) plane. The integration over k is performed with the variables ?/ = p2 ■ 1, 
X = k l, and the azimuthal angle (f)k. They are ready to be performed numerically. 

In contrast with Ref. |l] where the integrals contained in the parts corresponding to Eqs. (fl?)) 
and (fT5l) -namely, Eq. (35) of this reference- were most of them new and required a substantial 
effort to be performed analytically, all the integrals contained in Eqs. (fT4l) and (fT5l) have already 



been performed analytically in our previous work. It^s here where the parallelism between the 
present analysis and the one for the Si • I calculation [2] becomes very useful. One can compare 
Eqs. (fT4l) and ([T?] ) with their counterparts Eqs. (34) and (35) of this last reference and observe that 
the integrals with the common factor correspond to the integrals in our p^, while those with 
the common factor D3 correspond to our p^. Looking at the analytical results of the integrals in 
those Eqs. (34) and (35), given in Eqs. (46) and (47) of this same reference, we can establish the 
following connection 



and 



P3 



I \Sl-l 



(P2 + pT 



(18) 



=ip[+pir\ (19) 

where p[, . . . , p\ are given explicitly in Eqs. (48) and (49) of such Ref. |20. Here the upper labels 
■ 1 and Si • 1 are introduced to stress this correspondence. Using these results, and after some 
rearrangement to get somewhat more compact expressions, the analytical forms of Eqs. (fT4l) and 
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(fTSl) are 



P2 

P3 = y 



E\Y^ - F3) - le^E + + im2[2(l - - 5^3] + ^(3i?' - 2/^)^4 



I 



-^Eie, - (1 - 0')^e, + ^E^^7 + V + ^'^10 - o^^i4 - ^(4^ + ^")^o + ^ 



P2 



4 



2E 



,(20) 



(21) 



The explicit forms of the functions 6'j, 70, ?7o, Cij^ and Zi need not be reproduced here. They are 
all found in Ref. 

Collecting the virtual and bremsstrahlung RC our final result is 



dT (A- B^e-V) =dn{ U[) + -6/ 



TT 



sf . i 



(22) 



where the explicit forms of 0/ and A'q coincide with $1 of Eq. (54) and A'^ of Eq. (Bl) of Ref. ||2I], 
respectively. Our new results are A\^^ of Eq. di) and 



Qu = B'[(t> + Io{E, E2)] + B"<P' + D,ps + 



(23) 



where all the entries are defined above. The final result is infrared convergent. 

We should recall that the practical application in the Monte Carlo analysis may be to use the 
RC in the form 



Qi = atfl + hifigi + ctgl, 



(24) 



and to calculate the numerical values of the coefficients a^, hi, Ci throughout the Dalitz plot in 
the form of arrays. Such arrays would be fed into the Monte Carlo simulation as a matrix 
multiplication. This procedure should save a substantial computer effort. In Eq. ((24)) the index 
is i = /, //. 

The present results complement the ones of Ref. yj]. The RC to the Dalitz plot when the sf • 1 
angular correlation is observed have been obtained both in the rest frame of the emitted baryon B 
and in the rest frame of the decaying baryon A. They cover the three-body region of the Dalitz plot, 
they can be used for all charge assignments of A and B IVl], and / may be e^, /x^, or r^. They are 
presented in a form which is not compromised to fixing the values of the form factors at prescribed 
values. They provide a good approximation to RC of medium- (several tens of thousands) and low- 
(several thousands) statistics experiments in light- and heavy-quark baryon semileptonic decays, 
respectively. As a final remark, let us stress that even if we use the same notation of our previous 
work, the expressions here apply only to the present case and there should arise no confusion. 
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